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The Navier-Stokes order hydrodynamic equations for a low density granular mixture obtained 
previously from the Chapman-Enskog solution to the Boltzmann equation are considered further. 
The six transport coefficients associated with mass and heat flux in a binary mixture are given as 
functions of the mass ratio, size ratio, composition, and coefficients of restitution. Their quantitative 
variation across this parameter set is demonstrated using low order Sonine polynomial approxima- 
tions to solve the exact integral equations. The results are also used to quantify the violation of 
the Onsager reciprocal relations for a granular mixture. Finally, the stability of the homogeneous 
cooling state is discussed. 

I. INTRODUCTION 

The relevance and context of a hydrodynamic description for granular gases remains a controversial topic. At 
sufficiently low density the origin of hydrodynamics can be studied from Boltzmann kinetic theoryi^ To obtain 
hydrodynamics from the Boltzmann kinetic equation the essential assumption is the existence of a "normal" solution, 
defined to be one for which all space and time dependence occurs through the macroscopic hydrodynamic fields. 3 This 
solution, together with the macroscopic balance equations, leads to a closed set of hydrodynamic equations for these 
fields. The Chapman-Enskog method provides a constructive means to obtain an approximation to such a solution 
for states whose spatial gradients are not too large, tn this general context, the study of hydrodynamics for granular 
gases is the same as that for normal gases. 

The details of the Chapman-Enskog method have been carried out for a one granular component gas to obtain 
the Navier-Stokes order hydrodynamic equations, together with exact integral equations determining the transport 
coefficients occurring in these equations. 4 " These integral equations have been solved approximately using Sonine 
polynomial expansions and compared to both Direct Monte Carlo simulation of the Boltzmann equation and molecular 
dynamics simulation of the gas£ Good agreement is obtained even for relatively strong degrees of dissipation. The 
results support the formal theoretical analysis and the claim that hydrodynamics is not limited to the quasi-elastic 
limit. 

The analysis for multicomponent granular gases is much more complicated than for a one component gas. Most of 
the previous attempts^ were made for nearly elastic spheres where the equipartition of energy can be considered as 
an acceptable assumption, tn addition, according to this level of approximation, the inelasticity is only accounted for 
by the presence of a sink term in the energy balance equation, so that the expressions for the transport coefficients 
are the same as those obtained for normal fluids. 7 " However, the failure of energy equipartition in multicomponent 
granular gases- has been also confirmed by computer simulations^ and even observed in real experiments. 10 Although 
the possibility of nonequipartition was already pointed out many years ago^i it has not been until recently that a 
systematic study of the effect of nonequipartition on transport has been carcfuly analyzed. In this context, Garzo and 
Duftj* 1 ^ have carried out a derivation of the Navier-Stokes hydrodynamic equations for a binary mixture at low-density 
that accounts for nonequipartition of granular energy. These equations and associated transport coefficients provide 
a somewhat more stringent test of the analysis since the parameter space is much larger. There are now many more 
transport coefficients, given as functions of the three independent coefficients of restitution, size ratio, mass ratio, and 
composition. As in the one component case, explicit expressions for the transport coefficients requires also to consider 
Sonine polynomial expansions. The accuracy of this approach has been confirmed by comparison with Monte Carlo 
simulations of the Boltzmann equation in the cases of the shear viscosity^ 3 - and the tracer diffusion^ coefficients. 
Exceptions to this agreement are extreme mass or size ratios and strong dissipation, although these discrepancies 
between theory and simulation diminish as one considers more terms in the Sonine polynomial approximation^ 

Since the dependence of the shear viscosity coefficient on the parameters of the mixture (masses, sizes, concentration, 
coefficients of restitution) has been widely studied in a previous work^ a primary objective here is to demonstrate 
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the variation of the six transport coefficients associated with the mass and heat flux in this parameter space, using 
the same Sonine polynomial approximation as was found applicable for the one component gas. To set the context 
for these quantitative results two qualitative and potentially confusing issues are briefly noted at the outset. 

There is some ambiguity regarding the hydrodynamic temperature in a mixture since temperatures for each species 
can be defined in addition to the global temperature. What should be the hydrodynamic fields? For normal gases 
the answer is clear. These are set by the slow variables (at large space and time scales) associated with conserved 
quantities. In addition to species number and momentum, only the total kinetic energy is conserved, so only one 
global temperature occurs as a hydrodynamic field. The energy is no longer conserved for granular gases, but it 
remains a slow variable if the cooling rate is not too large. Thus, in this case as well, only the global temperature 
should appear among the hydrodynamic fields. 

Nevertheless, the species temperatures play a new and interesting secondary role. For a normal gas, there is a rapid 
velocity relaxation in each fluid cell to a local equilibrium state on the time scale of a few collisions. Subsequently, 
the equilibration among cells occurs via the hydrodynamic equations. In each cell the species velocity distributions 
are characterized by the species temperatures. These are approximately the same due to equipartition, and the 
hydrodynamic relaxation occurs for the single common temperature^ A similar rapid velocity relaxation occurs for 
granular gases in each small cell, but to a universal state different from local equilibrium and one for which equipartition 
no longer occurs. Hence, the species temperatures Tj are different from each other and from the overall temperature 
T of the celli& Nevertheless, the time dependence of all temperatures is the same in this and subsequent states, 
Ti(t) = jjT(t). This implies that the species temperatures do not provide any new dynamical degree of freedom. 
They still characterize the shape of the partial velocity distributions and affect the quantitative averages calculated 
with these distributions. The transport coefficients for granular mixtures therefore have new quantitative effects 
arising from the time independent temperature ratio 7$ for each species* 12 This dependence is illustrated here as well. 

In some earlier works^i*i^ additional equations for each species temperature have been included among the hydro- 
dynamic set. This is possible since the overall temperature is determined from these by T(t) = ^2iXiTi(t), where Xi 
is the concentration of species i. However, this is an unnecessary complication, describing additional kinetics beyond 
hydrodynamics that is relevant only on the time scale of a few collisions. As described above (and supported by 
molecular dynamics simulations^) the dynamics of all temperatures quickly reduces to that of T(t). The remaining 
time independent determination of the ji follows directly from the condition that the cooling rates of all temperatures 
be the same. 

A second confusing issue is the context of the Navier-Stokes equations considered here. The derivation of the Navier- 
Stokes order transport coefficients does not limit their application to weak inelasticity. For this reason the results 
reported below include a domain of both weak and strong inelasticity, 0.5 < a < 1. The Navier-Stokes hydrodynamic 
equations themselves may or may not be limited with respect to inelasticity, depending on the particular states 
considered. The derivation of these equations by the Chapman-Enskog method assumes that relative changes in 
the hydrodynamic fields over distances of the order of a mean free path are small. For normal gases this can be 
controlled by the initial or boundary conditions. It is more complicated for granular gases. In some cases (e.g., steady 
states such as the simple shear flow problem^) the boundary conditions imply a relationship between the coefficient 
of restitution and some hydrodynamic gradient - the two cannot be chosen independently. Consequently, there are 
examples for which the Navier-Stokes approximation is never valid or is restricted to the quasi-elastic limitiii However, 
the transport coefficients characterizing the Navier-Stokes order hydrodynamic equations are well-defined functions 
of a, regardless of the applicability of those equations (e.g., the coefficients of a Taylor series for some function f(a, x) 
in powers of x may be defined for all a at each order, even when the series cannot be truncated at that order). 

The plan of the paper is as follows. First, in Sec. [H] the hydrodynamic equations and associated fluxes to Navier- 
Stokes order are recalled, and the expressions for the transport coefficients for heat and mass transport are given 
to leading Sonine polynomial approximation. The elastic limit is discussed to aid in the interpretation of these 
expressions. Next, in Sec. lIIII the results for these six coefficients are illustrated for a common coefficient of restitution 
and same size ratio as functions of a at a composition x\ = 0.2 for several values of the mass ratio. With the exception 
of thermal conductivity, the deviations from normal gas values are largest at small a and large mass ratio. The usual 
Onsager relations among these coefficients for normal gases is then noted and tested for the granular gas in Sec. 
IIVI Since the underlying basis for these relations (time reversal symmetry) no longer holds for granular systems, the 
expected violation is demonstrated as a function of a for the same conditions. The stability of a special homogeneous 
solution to the mixture hydrodynamic equations is then studied, and some comments on the implications of the 
instability found are offered. Finally, the results are summarized and discussed in the last section. 
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II. BOLTZMANN KINETIC THEORY FOR THE MASS AND HEAT FLUXES 

We consider a binary mixture of inelastic, smooth, hard spheres of masses m\ and 77i 2 , and diameters <j\ and 02 • 
The inelasticity of collisions among all pairs is characterized by three independent constant coefficients of normal 
restitution a n , a22, and a 12 = a^x, where a.y is the coefficient of restitution for collisions between particles of species 
i and j. It is assumed that the density of each species is sufficiently low that their velocity distribution functions are 
accurately described by the coupled set of inelastic Boltzmann kinetic equations. The precise form of these equations 
is given in Ref. 12 and will not be required here. These equations imply the exact macroscopic balance equations for 
the particle number density of each species, (r, t), flow velocity u (r, t), and temperature, T (r, t)^ 

Dtni + ntf-ul — =0 , i = 1, 2 (1) 

mi 

D t u + p _1 VP = 0, (2) 

OtT--E — + |-(V-q+P:Vu) = -CT. (3) 

i 

In the above equations, D t = dt + u • V is the material derivative, p = mini + mini is the total mass density, ji 
is the particle number flux for species i, q is the heat flux, P is the pressure tensor, and £ is the cooling rate. For 
the two component mixture considered here there are six independent fields, n±, n%, T, u. To obtain a closed set 
of hydrodynamic equations, expressions for j^, q, P, and C must be given in terms of these fields. Such expressions 
are called "constitutive equations" . It is convenient to give these constitutive equations in terms of a different set 
of experimentally more accessible fields, x\, p, T, u, where x\ = n\/ (n± + ni) is the composition of species 1, and 
p = (ni + n 2 ) T is the hydrostatic pressure. This is simply a change of variables, so that Eqs. CJ-JSJ become 

D tXl + ^ V-ji=0, (4) 

2 

D t p + pW ■ u + - (V • q + P : Vu) = -(p, (5) 

D t u+p- 1 VP = 0, (6) 

AT--£-^ + — V-q+P:Vu ) = -(T. (7) 
n i! — ' mi 3n 

i 

The constitutive equations up to the Navier-Stokes order have been obtained from the Boltzmann equation in Ref. 
H2I with the results 

ji = - ( -^-JDVxx-^DpVp-^D'VT, j 2 =-ji, (8) 
q = -T 2 D"\7 Xl - LVp - AVT, (9) 
Pkt = p5ki - V ^V f ufe + V k u e - ^S ke W ■ u^j , (10) 

C = Co + G(V 2 ) (11) 

The transport coefficients {D, D p , D', D" , L, A, 77} verify a set of coupled linear integral equations which can be solved 
approximately by using the leading terms in a Sonine polynomial expansion. This solution provides explicit expressions 
for the transport coefficients in terms of the coefficients of restitution and the parameters of the mixture (masses, 
sizes, and composition). The above expressions for mass and heat fluxes can be defined in a variety of equivalent 
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ways depending on the choice of driving forces used. For systems with elastic collisions, the specific set of gradients 
contributing to each flux is restricted by fluid symmetry, Onsager's relations (time reversal invariance), and the form 
of entropy production*^ In this case, one usual representation leads to the mass and heat fluxes proportional to 
(V/Xj)x and VT, where \ii (defined below) is the chemical potential per unit mass. However, for inelastic systems 
only fluid symmetry holds and so there is more flexibility in representing the fluxes and identifying the corresponding 
transport coefficients. In particular, a third contribution proportional to Vp appears in both fluxes. Some care is 
required in comparing transport coefficients in different representations using different independent gradients for the 
driving forces. 

The cooling rate to lowest order in the gradients is C = Co( x i,P, T). There are no contributions to first order in the 
gradient for the low density Boltzmann equation. The general form including second order gradient contributions is 
displayed in Appendix A. These second order terms have been calculated for a one component gas^ and found to be 
very small. Here, these second order contributions to the cooling rate will be neglected. 

Substitution of the Navier-Stokes constitutive equations, (JSJ-lJTSJ, into the exact balance equations, (0}-(0), gives 
the Navier-Stokes hydrodynamic equations for a binary mixture 



D t xi 



V • 



V P 



DVx! + -D p S7p+ ^D'VT) , 



T 



(12) 



■V ■ (T 2 D"X7 Xl + LVp + AVT) 

t 

2 / 2 



(13) 



(D t + C)T+-pV-u 



T mi — ra\ 
n TO1TO2 



V ■ 



DVxi + -D p Vp + -D'VT 
P V T 



+— V ■ (r 2 D"Vxi + LVp + AVT) 

2 / 2 \ 

+ 3 -r ? yJi u k + V fc Uf - -8 ke V ■ u I V e u k , 



(14) 



D t u e + p ^ep = p ( V t u k + V k u t - -5 ke V ■ u 



(15) 



For the chosen set of fields n = p/T and p = p [{mi — 777.2 ) %i + 1^2] /T. These equations are exact to second order in 
the spatial gradients for a low density Boltzmann gas. 

A. Mass flux 

The mass flux contains three transport coefficients, the diffusion coefficient D, the pressure diffusion coefficient D pi 
and the thermal diffusion coefficient D' . Explicit expressions for these were obtained in Ref. using a first Sonine 
approximation. Dimensionless forms are defined by 



D 



PT 



777 1 7772 ^0 



-D* 



D - -D* 

P ~ pvo p ' 



D> = ^D> 
pvo 



(16) 



Here, uq = ^/nna^vo, o\2 — (<7i + ^2)/^, and vq — y / 2T(mi + 7772) / '77717712 is a thermal velocity defined in terms of 
the temperature T of the mixture. The explicit forms are then 



D* 



d 
dx\ 



£171 



9C 
dx\ 



— ! n 

2v* 1 p 



(17) 



d; = Xi 



71 



m(i + s) 



l + [x5 



* 3 A * C 2 

v* C + - — 

2 s 2v* 



(18) 
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D'* = -^—D* (19) 
2v* p 



In these equations, 



T\ 1 T 2 1 , . 

71 = 7^7 = 7— 7 7T i 72 = 777 = "j— 7 7T; (. 2U ) 

T 1 + 2:1(7— 1) T 1+0:1(7- 1) 



where /i = mi/ni2 is the mass ratio, fuj — m.i/(m,i + mj), 5 — Xi/x%, and 7 = T1/T2. The detailed forms for the 
temperature ratio 7, dimensionless collision rate v* , and dimcnsionlcss cooling rate £* are given in Appendix A. Since 
ji = — J2 and Vcci = — VX2, it is expected that D* should be symmetric with respect to interchange of particles 1 
and 2 while D* and D'* should be antisymmetric. This can be easily verified by noting that £171 + £272 = 1- 
In the case of elastic collisions ctij = 1, £* = 0, 7 = 1, and Eqs. (|17|l 1191) become 

D * = l, F7 ' d;^ Xi ^—4d*, D'* = 0. (21) 

81-^12(1-5) p 1 + pt<5 v ; 

These coincide with known results obtained for elastic collisions in the first Sonine approximation^ Recently, it has 
been shown that the estimate given by the first Sonine approximation for the diffusion coefficient D (in the very 
dilute concentration limit X\ — > 0) compares quite well with Monte Carlo simulations, 14 : except for the cases in which 
the gas particles are much heavier and/or much larger than impurities. For these extreme cases, the second Sonine 
approximation to D improves the accuracy of the kinetic theory results^ 



B. Heat flux 



The heat flux requires going up to the second Sonine approximation. The transport coefficients D" , L, and A 
appearing in the heat flux 10 are given bj^ 



D" 



2 (mi + m 2 )v 



M12 



M21 
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71 
M12 



72 
M21 



D* 



(22) 



L = — - 



2 (mi + m 2 )^o 



M12 



M21 



f 7i 72 
VM12 M21 



(23) 



A = --- 



nT 



2 (mi + m 2 )^o 



^i 3 Ai + ^7| A2 



M12 



M21 



7i 
M12 



72 

M21 



(24) 



where the expressions for the (dimensionless) Sonine coefficients {d",£i,Ai} are displayed in Appendix lAl In Eqs. 
(|2*2*)l - l|2*4l it is understood that the coefficients D*, D* and D'* are given by Eqs. (JTTJhJTHJ, respectively (first Sonine 
approximation). As expected, our results show that D" is antisymmetric with respect to the change 1 <-> 2 while L 
and A are symmetric. Consequently, in the case of mechanically equivalent particles (mi = m.2 = m, o\ = a% = a, 
onj = a), the coefficient D" vanishes. 

An equivalent representation is given in terms of the heat flow 3 q defined as 



5„^-^ i, 5 „ mo — nil , 

2 z — ' rrii 2 m\mi 



(25) 



where in the second equality use has been made of the requirement ji = — J2. The difference between q and 3 q is a 
heat flow due to diffusion. In addition, for elastic collisions, 3 q is the flux conjugate to the temperature gradient in 
the form of the entropy production where the contribution coming from the mass flux couples only to the gradient 
of the chemical potentials. The thermal conductivity in a mixture is generally measured in the absence of diffusion, 
i.e., when ji = 0. To identify this coefficient, we have to express J q in terms of ji, VT, and Vp. The corresponding 
coefficient of VT defines the thermal conductivity 18 According to Eq. |JSJ, the gradient of mole fraction Vx\ is 



Vzi 



P 



D' 



7J1 



m\m2TiD m\m2P D 



-VT 



mim2np D 



(26) 
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The expression of J g is obtained by substituting Eq. into Eq. (|25|) and eliminating Vx± by using the identity i|26|l . 
Thus, the heat flow is given by 



where 



K = A 



pp 



miTO2?T. 



r K T ji - L p \7p, 



p 2 T 2 D ,t D , 

m\ui2y D 



(27) 



(28) 



T D" 5n 
n T = T— - --{m 2 - mi), 
D 2 p 



(29) 



p 2 T D"D p 
n 2 mirri2 D 



(30) 



As in the elastic case, the coefficient k is the thermal conductivity while kt is called the thermal-diffusion factor or 
Dufour coefficient. There is a new contribution proportional to Vp not present in the elastic case that defines the 
transport coefficient L p . 

For elastic collisions, L p = Oi^ and the expressions derived here for k and kt coincide with those obtained for a gas- 
mixture of elastic hard spheres. 3 Furthermore, in the case of mechanically equivalent particles, the Dufour coefficient 
kt vanishes as expected and the heat flow l|27|) can be written as 



-kVT — /iVn, 



(31) 



where 



k = k + nL 



25 / 



mT 



P ~ 32 V 7T 



1/2 



, 1 + c 



(32) 



75 T (mT\ 
32 n \~J 



1/2 



(33) 



r(i 



, 33, , 19 -3a 

1 -I (f -a)-\ c 

f6 v ; 1024 



(34) 



32(f-a)(f-2a 2 ) 
81 - f7a + 30a 2 (f - a)' 



12 



[I -a* 



32 



(35) 



Note that in writing Eq. I|31|l use has been made of the relation Vp = riVT + TWn. Equations (|31|I - H35|I are the 
same as those obtained for the one component granular gas^ This confirms the relevant known limiting cases for the 
granular mixture results described here. 



III. TRANSPORT COEFFICIENTS 



The transport coefficients depend on many parameters: {x\, T, m\/m-z, ci/o^, an, a22, oiyi}- This complexity exists 
in the elastic limit as well, so the primary new feature is the dependence on the coefficients of restitution ctij being 
different from unity. To illustrate the differences between granular and normal gases the transport coefficients are 
normalized to their values in the elastic limit. Then, the dependence on the overall temperature scales out. Also, 
only the simplest case of a common coefficient of restitution (an = a22 = ct\2 = ct) and common size u> = a\/<72 = 1 
is considered. This reduces the parameter set to three quantities: {mi/m2, x±, a}. 

In Figs. Hill we plot the above transport coefficients as functions of the coefficient of restitution a for x\ = 0.2, 
and several values of the mass ratio /i. It is understood that all coefficients have been reduced with respect to their 
elastic values, except in the cases of D' and L p since both coefficients vanish for elastic collisions. In these latter two 



7 



t 1 1 1 1 1 r 




I i i i i i i i i i I 

0.5 0.6 0.7 0.8 0.9 1.0 

a 



FIG. 1: Plot of the reduced mutual diffusion coefficient D(a)/ D(l) as a function of the coefficient of restitution a for x\ = 0.2, 
uj — 1 and /i = 0.5 (a), fi = 1 (b), and /i = 4 (c). The dashed lines correspond to the approximation ci = C2 = 0. 
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0.5 0.6 0.7 0.8 0.9 1.0 
a 

FIG. 2: Plot of the reduced pressure diffusion coefficient D p (a)/D p (l) as a function of the coefficient of restitution a for 
xi = 0.2, uj = 1 and fj, = 0.5 (a), /j, = 2 (b), and fi — 4 (c). 

cases, we have considered the reduced coefficients D'* defined by Eq. (|17fl and L* — ~ [(5/2)TVo/(toi + m^)] -1 L p . 
Figure n shows the mutual diffusion coefficient as a function of a for three mass ratios fj, = 0.5, 1, and 4. There is a 
monotonic increase of the coefficient with decreasing a in all cases. Moreover, that effect increases as the mass of the 
dilute species increases. This is consistent with an observed singular behavior in the extreme case of tracer diffusion 
for a massive particle^ The velocity distributions in a granular gas are no longer Maxwellian^ and the difference is 
measured by the coefficients a that appear in the expressions for the transport coefficients^ The dashed curves in 
Fig. ^ correspond to c\ = — > 0, the Maxwell limit. In this case it is seen that the distortion of the Maxwellian is 
not very important for this transport coefficient (see however, discussion of Fig. 0] below) . Figure [21 shows that the 
pressure diffusion coefficient has a very similar behavior. The thermal diffusion coefficient vanishes in the elastic limit 
and remains small and slightly negative when the dilute species has small mass ratio, as illustrated in Fig. EI However 
as the mass ratio becomes large it becomes large and positive for strong dissipation. The effect of different species 
temperatures is also shown on this graph. The dashed curves correspond to setting 7 = X1/T2 — + 1. This is seen to 
yield large errors, particularly in the case of large mass ratio (temperature differences are greater for mechanically 
different particles), indicating the real quantitative effect of two different species temperatures in granular gases. 

The thermal conductivity is shown in Fig. with the same monotonic increase with increasing dissipation. The 
dashed lines (cj = 0) indicate a significant effect of the distortion of the reference distribution function from its 
Maxwellian form. Presumably, this is due to the fact that the thermal conductivity k depends on a higher velocity 
moment than the mutual diffusion coefficient D and is more sensitive to the larger distortions at higher velocities. 
We also observe that there is little mass dependence when the dilute species is lighter. However, when the dilute 



8 



0.5 
0.4 
0.3 
D 0.2 
0.1 
0.0 
-0.1 

0.5 0.6 0.7 0.8 0.9 1.0 
a 

FIG. 3: Plot of the reduced thermal diffusion coefficient D *(a) as a function of the coefficient of restitution a for x\ — 0.2, 
u = 1 and n = 0.5 (a), \jl = 2 (b), and [s, = 4 (c). The dashed lines correspond to 7 = f . 



2.00 
1.75 

f — s 

'a 

V 1.25 
1.00 



0.5 0.6 0.7 0.8 0.9 1.0 
a 

FIG. 4: Plot of the reduced thermal conductivity coefficient k(o)/k(1) as a function of the coefficient of restitution a for 
x\ = 0.2, u> = 1 and = 0.5 (a), pi = f (b), and fi — 4 (c). The dashed lines correspond to the approximation ci = C2 = 0. 

species is more massive there is a significant decrease in the thermal conductivity, opposite to the case of diffusion. 
In contrast, the Dufour coefficient does have a dependence on the mass ratio more like diffusion (Fig. 0. Finally, the 
coefficient L* is illustrated in Figure El Again there is weak dependence on the mass ratio until it becomes larger for 
the dilute species. 

In summary, the mass and heat flux transport coefficients for a granular mixture differ significantly from those for a 
normal gas mixture even at moderate dissipation. In most cases (thermal conductivity is an exception) the differences 
increase with decreasing a, depend weakly on the mass ratio when the dilute species (x\/x2 < 1) is lighter than the 
excess species (mi/ma < 1) but increase significantly in the opposite case (mi/ms > 1). 

IV. ONSAGER'S RECIPROCAL RELATIONS 

In the usual language of the linear irreversible thermodynamics for ordinary fluids^ the constitutive equations for 
the mass flux © and heat flow l|27|l is written 

U = E L « (ir) t L ^YT - C ^P, (36) 





J i I i I i L 



0.5 0.6 0.7 0.8 0.9 1.0 
a 



FIG. 5: Plot of the reduced Dufour coefficient kt(q)/kt(1) as a function of the coefficient of restitution a for x\ = 0.2, lo — 1 
and = 0.5 (a), /i = 2 (b), and fi = 4 (c). 




a 



FIG. 6: Plot of the reduced coefficient Lp(a) as a function of the coefficient of restitution a for X\ = 0.2, u = I and fi — 0.5 
(a), jj, — 1 (b), and fi = 4 (c). 



J q = -L qq VT - L qi ( ^ ) _ - C;Vp, (37) 



V/ ''^ -— Vln(*iP)» (38) 



where 

H- , _ 

T J T mi 

\Xi being the chemical potential per unit mass. Here, the coefficients are the so-called Onsager phenomenological 
coefficients. For normal fluids, Onsager showed^ that time reversal invariance of the underlying microscopic equations 
of motion implies important restrictions on the above set of transport coefficients 

Lij = Lji 1 Liq = Lqi, C'p — Cp = 0. (39) 

The first two symmetries are called reciprocal relations as they relate transport coefficients for different processes. 
The last two are statements that the pressure gradient does not appear in any of the fluxes even though it is admitted 
by symmetry. Even for a one component fluid, Onsager's theorem is significant as it leads to Fourier's law for the 
heat flow rather than IpU , i.e. ~p = 0. Since there is no time reversal symmetry for granular fluids, Eqs. (|39|) cannot 
be expected to apply. However, since explicit expressions for all transport coefficients are at hand, the quantitative 
extent of the violation can be explored. 
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To make connection with the previous sections it is first necessary to transform Eqs. I|36|l - (|38(l to the variables x\, 
p, T. Since Vx x = -Vz 2 , Eq. (J3BJ) implies 



(V/^i)t - (V/i 2 ) T 



np 
P1P2 



n\n 2 

vx\ H (m 2 — mijvlnp 

np 



The coefficients Ly then can be easily obtained in terms of those of the previous sections. The result is 



(40) 



in 
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= 5 D, Li q = pTD , 



(41) 
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T 2 pip 2 „„ 5 Tpip 2 



-D 



np 
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D 



(42) 



G p = -L>„ H™2 - "ii ) A 

p pp 



(43) 



C' =L - 5 T m2 ~ mi C - ni " 2 ^ L> 
2 p vn\Tfi 2 p npp 



T 2 (m 2 — mi)D" . 



(44) 



The Onsager's relation Li 2 = L 2 i is already evident since the diffusion coefficient D is symmetric under the change 
1 <-* 2 , as discussed following Eq. I|20|l 

Imposing Onsager's relation Li g = L g i yields 



„/, 5 n , v no „. 

2 rp t pip2 

while the condition C p — CL — leads to the following additional requirements 



(45) 



U„ = — ^-(m 2 - mi)D, 



(46) 



5T 

--(771! - m 2 

2 P 



(m 2 - 77ii)D 



m\m 2 



"A; 



71,0 



(47) 



Since the relations l|45(l - (|47|l involve transport coefficients that have been determined in the first Sonine approximation, 
we restrict here our discussion to this level of approximation. In this case, d'[ = £ j = A, = so that Onsager's theorem, 
Eqs. iP^jl - PTIl . gives the conditions 



(48) 



Dl 



Xl TTa — x\ D =0 ' 



(49) 



R(aij) 



1 + M 



(71 -M72) Q(Qfy) = 0. 



(50) 



In the elastic limit, the reduced coefficients D* and D* are given by Eq. (|19|l and these conditions are verified. 
Also, for mechanically equivalent particles with arbitrary a, 7, = 1 and D* = so that P(a) — Q(a) = R(a) = 0. 
Nevertheless, beyond these limit cases, Onsager's relations do not apply (as expected). At this macroscopic level 
the origin of this failure is due to the cooling of the reference state as well as the occurrence of different kinetic 
temperatures for both species. Figures El and show the dependence of the quantities P, Q, and R, respectively, 
on the (common) coefficient of restitution = a for mass ratios p = 0.5, 2, and 4. 




FIG. 9: Plot of R(a) as a function of a for xi = 0.2, uj = 1 and /i = 0.5 (a), fj, = 2 (b), and /i = 4 (c). 
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V. LINEARIZED HYDRODYNAMIC EQUATIONS AND STABILITY 

In contrast to normal fluids, the Navicr-Stokes hydrodynamic equations l|12l) - (|15[) have non-trivial solutions even 
for spatially homogeneous states, 

d t x 1H = = d t u Ht , (51) 

[dt + C(xih,Th,Ph)}T h = 0, [dt + {(x 1H ,T H ,PH)}PH = 0. (52) 

where the subscript H denotes the homogeneous state. Since the dependence of the cooling rate £ {x\h,Th,Ph) on 
xiHi Th,Ph is known (see Appendix A), these first order nonlinear equations can be solved for the time dependence 
of the homogeneous state. The result is the familiar Haff's law cooling law for T(t) at constant density. 2 . As 
discussed above, each species temperature also has the same time dependence but each with a different valued In 
this section, the hydrodynamics for small initial spatial perturbations of this homogeneous cooling state (HCS) is 
discussed. For normal fluids such perturbations decay in time according to the hydrodynamic modes of diffusion 
(shear, thermal, mass) and damped sound propagation. The analysis is for fixed coefficients of restitution different 
from unity in the long wavelength limit. It will be seen here that the corresponding modes for a granular fluid are 
then quite different from those for a normal fluid. An alternative analysis with fixed long wavelength and coefficient 
of restitution approaching unity leads to usual normal fluid modes. Thus, the nature of hydrodynamic modes is non 
uniform with respect to the inelasticity and the wavelength of the perturbation. 

Let 5y a (r,t) = y a (r,t) — yHaif) denote the deviation of {xi,u,T,p} from their values in the HCS. If the initial 
spatial perturbation is sufficiently small, then for some initial time interval these deviations will remain small and 
the hydrodynamic equations (|12|> -H15 [) can be linearized with respect to 5y a (r,t). This leads to partial differential 
equations with coefficients that are independent of space but which depend on time since the HCS is cooling. As in 
the one component case^Si this time dependence can be eliminated through a change in the time and space variables, 
and a scaling of the hydrodynamic fields. We introduce the following dimensionless space and time variables: 

rt 

dt'v H {t'), r'=r/l H , (53) 



where vqh = \J1Tn(m\ + m2)jm\mi is the thermal velocity introduced above, In = 1/ is an effective mean 
free path, and vh{^) — v oh(^)/^H is the effective collision frequency. The dimensionless time scale is therefore an 
average number of collisions up to the time t. A set of Fourier transformed dimensionless variables are then defined 
by 

<Wr) = J dr' e-^ r 'Sy a (r',T), (54) 

<faik(r) <5u k (r) . . 

Pk r = , w k (r) = r -, (55) 

XlH Vqh(T) 

Ok t = ^TTT i n k(T = rr- (56) 

In terms of these variables the linearized hydrodynamic equations for the set {pk, wk, #k, Hk} separate into a degen- 
erate pair of equations for the two transverse velocity components Wk± (orthogonal to k) 

'' " 1 V*k 2 ) Wki = 0, (57) 



(M$ + ikM^ + k 2 M$) 5z k p(r), (58) 



^dr 2 

and a coupled set of equations for /9k,#k,Hk, and the longitudinal velocity component Wk| (parallel to k) 

dSzk a (T 

dr 

where now Sz^ a (T) denotes the four variables (pk, ^k 5 Hk, Wuii) ■ The matrices in this equation are 

/ \ 



£ J 



(59) 
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M« = 











\ 

_ 2 




















~ 3 


V° 





1 /*12 

2 Xi[l-\-X2 


/ 



(60) 



V 



1 fJ,X 1 +X 2 TJ* 

2 1 + 
2 rV'* 1 — ^ rj* 

2(1+/*) 





__ 1 fix 1 +x 2 jjl* 
2xi 1+/* 

1 — rv* _2 \* 

2(1+/*)^ 3 A 
2 \ * 
~3 A 





1 ^a;i+x 2 7-i* 

2a;i 1+jU P 
2 r * i 1 — /* n* 
3"^ + 2(1+/*) "^P 
2 r* 
3^ 








(61) 



In these equations, Xi = rim/nHi C* = Ch/^h and the reduced transport coefficients D* , D*, and D'* are given by 
Eqs. (|17fl - l|19l 1. respectively. Moreover, the reduced Navier-Stokes transport coefficients are 



V = 



(62) 



_ v H T H D" 



D = 



n H VQH 



L* = 



v H L 



J 0H 



A* = 



(63) 



where pn — minm + m 2 n 2H- 

It is instructive to consider first the solutions to these equations in the extreme long wavelength limit, k = 0. In 
this case, the eigenvalues or hydrodynamic modes are given by 



1 



0,0, 



2^ ' 2^ 



(64) 



where s n refers to the longitudinal modes. Two of the eigenvalues are positive, corresponding to growth of the initial 
perturbation in time. Thus, some of the solutions are unstable. The two zero eigenvalues represent marginal stability 
solutions, while the negative eigenvalue gives stable solutions. For general initial perturbations all modes are excited. 
These modes correspond to evolution of the fluid due to uniform perturbations of the HCS, i.e. a global change in 
the HCS parameters. The unstable modes are seen to arise from the initial perturbations WkJ_(0) or u>k||(0). The 
marginal modes correspond to changes in the composition at fixed pressure, density, and velocity, and to changes 
in Ilk — $k at constant composition and velocity. The decaying mode corresponds to changes in the temperature or 
pressure for Ilk = $k- The unstable modes may appear trivial since they are due entirely to the normalization of the 
fluid velocity by the time dependent thermal velocity. However, this normalization is required by the scaling of the 
entire set of equations to obtain time independent coefficients. 

At finite wave vectors, these instabilities give rise to real growth of spatial perturbations. The linear growth of the 
transverse modes is simply given by 



WkiW = w k +(0)exp(ic* - r\*k 2 )T. 
The instability for the two shear modes is removed at sufficiently large k > kj_, where 

1/2 



k c , = 



2ry* 



(65) 



(66) 



The wave vector dependence of the remaining four modes is more complex. This is illustrated in Fig. 1101 showing the 
real parts of the modes s (k) for = 0.9, u\ju2 = 1, x% = 0.2, and m%/m2 = 4. The k = values are those of l|64|l . 
corresponding to six hydrodynamic modes with two different degeneracies. The shear mode degeneracy remains at 
finite k but the other is removed at any finite k. At sufficiently large k a pair of real modes become equal and become 
a complex conjugate pair at all larger wave vectors, like a sound mode. The smaller of the unstable modes is that 
associated with the longitudinal velocity, which couples to the scalar hydrodynamic fields. It becomes negative at a 
wave vector smaller than that of Eq. (|66|l and gives the threshold for development of spatial instabilities. 

The results obtained here for the mixture show no new surprises relative to the earlier work for a one component 
gas^*2i with only the addition of the stable mass diffusion mode. Of course, the quantitative features can be quite 
different since there are additional degrees of freedom with the parameter set {x±, T, m\/m2, 01/021 . Also, the 
manner in which these linear instabilities are enhanced by the nonlinearities may be different from that for the one 
component case. 
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FIG. 10: Dispersion relations for a — 0.9, x\ = 0.2, u) = 1 and /j, = 4. 

VI. SUMMARY AND DISCUSSION 

The Navier-Stokes order hydro-dynamic equations have been discussed for a low density granular binary mixture. 
The form of the momentum flux is the same as for a one component gas, with only the value of the viscosity changed 
(see Appendix 0). Since the dependence of the viscosity on the parameters of the mixture has been widely explored 
in a previous paper^ attention has been focused here on the mass and heat fluxes and their associated transport 
coefficients. There is no phenomenology involved as the equations and the transport coefficients have been derived 
systematically from the inelastic Boltzmann equation by the Chapman-Enskog procedure. Consequently, there is 
no a priori limitation on the degree of inelasticity, size and mass ratios, or composition. For practical purposes, 
the integral equations determining the transport coefficients have been solved by truncated expansions in Soninc 
polynomials. This is expected to fail at extreme values of size or mass ratio^iil but the results are quite accurate 
otherwise. 

The hydrodynamic equations are the same as for a normal gas, except for a sink in the energy equation due 
to granular cooling, and additional transport coefficients in the mass and heat flux constitutive equations. The 
latter arise because the usual restrictions of irreversible thermodynamics no longer apply. These restrictions include 
Onsager reciprocal relations among various transport coefficients, and the extent to which these are violated has been 
demonstrated in Section 4. It has been verified that the results described here reduce to those for a normal mixture 
in the elastic limit, 3 and to those for a one component granular gas4 when the species are mechanically identical. 

As is the case for a normal gas, the hydrodynamic fields include only the global temperature even though two species 
temperatures can be defined. For a normal gas the species temperatures rapidly approach the global temperature 
due to equipartition. For the granular gas the species temperatures approach different values, but with the same time 
dependence as the global temperature. The transport coefficients have an additional dependence on the composition 
due to this time independent ratio of species temperatures. This has been illustrated in Fig.[3]for the thermal diffusion 
coefficient, where the effect is seen to be large for large mass ratio (implying very different species temperatures. 
Generally, it was seen that the deviation due to inelasticity is enhanced for greater mechanical differences between 
the species. 

The linear equations for small perturbations of the special homogeneous cooling solution to the hydrodynamic 
equations were obtained and discussed. In order to characterize the solutions in terms of modes, it is necessary to 
introduce dimensionless fields so that the time dependence of the reference HCS is eliminated. The resulting equations 
exhibit a long wavelength instability for three of the modes. This is quite similar to the case of a one component 
granular gas^Si and in fact the same modes are unstable here. The additional diffusion mode for two species behaves 
as for a normal fluid. The consequences of this instability for a binary mixture were not studied here. This entails 
an analysis of the dominant nonlinearities which has not been performed as yet. Since there are additional degrees of 
freedom would be interesting to see if the density clustering that occurs for a one component system is more complex 
here (e.g., species segregation). 
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APPENDIX A: SOME EXPLICIT EXPRESSIONS 



Navier-Stokes hydrodynamics retains terms up through second order in the gradients. As a scalar, the cooling rate 
has the most general form at this order given by 

C = Co + C«v • u + GV^i + CtV 2 T + C P V 2 p + Ctt (VT) 2 + ( xx (Van) 2 
+C PP (Vp) 2 + Ct x (VT) ■ (Van) + Ct p (VT) • (Vj») + ( px (Vp) ■ (Vz x ) 
+Cuu {ViUj) {ViUj) . (Al) 

For a low density gas the first order gradient term vanishes, 4 ( u = 0. However, for higher densities ( u is different 
from zeroi^ 2 - The second order terms have been left implicit in equations (|13|l and Ijl4|l for the temperature and the 
pressure. As noted in the text, these second order terms have been calculated for a one component fluid^ and found 
to be very small relative to corresponding terms from the fluxes. Consequently, they have been neglected in the 
linearized hydrodynamic equations (|58|l . 



1. Mass flux parameters 



The transport coefficients are expressed in terms of a number of dimensionless parameters. For completeness, they 
are listed here. In general they depend on the reference distribution functions in the Chapman-Enskog expansion, 
which are not Maxwellians. 8 The deviation of the reference distributions from their Maxwellian forms is measured by 



the cumulants Cj. As shown in Sec. IIIII while the influence of these coefficients on the transport coefficients is not 
quite important in the case of the coefficients associated with the mass flux (see Fig. Q f° r example), not happens 
the same in the case of the heat flux (see Fig. for example) where the influence of Cj is not negligible for strong 
dissipation. However, in order to offer a simplified theory the parameters given in this Appendix have neglected the 
corrections due to the cumulants ci. The full expressions for the transport coefficients can be found in Ref. la 
The temperature ratio 7 = T1/T2 is determined from the condition 



/-* /~* 

si S2 s 1 

where the dimensionless cooling rate, Q = Ci/^o { u o is the average frequency defined below (|lfc>|) ^ is 

cr = f^(^) 2 Mr 1/a (1-0 

+^X2iM,i (^r) ' i 1 + "12) e 2 1/2 [2 - M21 (1 + ai 2 ) (1 + 9)] . 



(A2) 



(A3) 



The expression for Q can be easily obtained by interchanging 1 <-» 2. Here, 9\ — 1/(^2171), #2 = 1/(^1272), 
Hij = mijirtii + rrij), S — x%/x2, and 6 — 61/82 = m/7- The temperature ratios ji are related to the temperature 
ratio 7 through Eqs. (|20|l . In the quasielastic limit (ojj -C 1), the temperature ratio 7 has the simple form 



1 



2^12^2: 



{{H12X1 - H21X2) (1 - CC12) 



1 

"71 



022 
C12 



"22) 



CTll 
Cl2 



- an) 



(A4) 



In this limit, the temperature ratio is a linear function of the hydrodynamic field x±. The dimensionless frequency v* 
appearing in the expressions of the transport coefficients associated with the mass flux is 



4 1 
0M21 



fj,S (\ 



1 + 5 



1/2 



)- 1/2 (l + a 12 ). 



(A5) 
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The cooling rate £* and frequency v* are also functions of the hydrodynamic field X\. However, all the parameters 
above are independent of the temperature and density. 



2. Heat and momentum flux parameters 

As in the last section effects due to the distortion of the reference Maxwellian are neglected. In the case of the heat 
flux q, Eq. ((?JJ, the transport coefficients D", L, and A are given by Eqs. (|22[) - l|24[l . respectively. By using matrix 
notation, the (dimcnsionless) Sonine coefficients d" , £i, and A^ verify the coupled set of six equations^ 



where X„i is the column matrix 



(A6) 



X 



Ai 
VA 2 / 



(A7) 



and A CT(T ' is the matrix 



A : 



^21 















-(SSL) 
\ dxi J P ; 



V\\ - 2C 

C/2 




-(SSL) 


C/2 



.(SSL) 

\ dxi J P : 





y\\ - C* 

V%1 



\ 

.(SSL) 

\ dxi Jp.T 



-C* 
^12 

-e J 



(A8) 



The column matrix Y n has the elements 



Y 1= D*[ na 



C* 



1 / 571 



Xi7f ) 7f \dx 



-D* T 2 



X27I 



1 / 972 



75 \dx 



(A9) 



la = D; ( ria 



1 

7i 



£>'* 712 - 



^i7i 



X271 
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F 6 = £>'* 7-21 - 

72 



Z27! 



(A10) 
(All) 



The dimensionless collision integrals T12, U11, and ^12 are given, respectively, by 



712 



4 / ^21 / _£! 



3V 2 Vcr 12 



-3/2., 2 \ 

7i (l-Oii) 



15 \ x\ 



(A12) 



■l5^ 2l(1 + ai2) VW+^ 



1 33 M 

1 + — (1-an) 

3/2 



(A13) 



2 M21 / 



3/2 



P2 



(A14) 
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In the above equations we have introduced the quantities^ 

A = 5{2[3 12 + 9 2 )+ f i 21 (9 1 +9 2 )[5(l-a 12 )-2(7a 12 -ll)p 12 9- 1 ] 

+ 18/3 2 2 9^ + 2/4 (2a 2 2 - 3a 12 + 4) 9^ 1 {9 1 + 9 2 f - 59 2 9^ 1 {9 1 + 9 2 ) (A15) 

B = b{2p 12 -e 1 )+ m (e 1 +9 2 )[5(l- ai2 ) + 2(7a 12 -n)f3 12 9 2 1 ] 

-IZPl^- 1 ~ 2/4 (2a 2 12 - 3a 12 + 4) 9 2 1 (9 1 + 6 2 ) 2 + 5(9 1 + 6 2 ) (A16) 

E = 2y 2 1 9^ 2 {9 l +9 2 f (2a 2 12 - 3a 12 + 4) (50! + 86 2 ) 

-IMnQBi + 2 ) [2/?i 2 6»r 2 (56»i + 80 2 )(7a 12 - 11) + 26 2 e^ 1 (29a 12 - 37) - 25(1 - a 12 )] 

+i8/3? 2 0r 2 (5^i + 89 2 ) + 2(i 12 e^ 1 {2he 1 + me 2 ) 

+59 2 9^ 1 (U9 1 + 69 2 ) - 5(0i + 9 2 )9, 2 9 2 {59 1 + 69 2 ) (A17) 

F = 2fi 2 21 9 2 2 {9 1 +9 2 ) 2 (2a 2 12 -3a 12 + 4) (8#i + 5# 2 ) 

-M2i(0i + 2 ) [2/3i 2 ^ 2 (8^ + 50 2 )(7a 12 - 11) - 29 1 9 2 1 (29a 12 - 37) + 25(1 - a 12 )] 
+ 18(3 2 12 9 2 2 {89 1 + 59 2 ) - 2/? 12 ^ 1 (660i + 250 2 ) 

+59 1 9 2 1 (69 1 + 119 2 ) - 5(6»x + 9 2 )9 2 1 (69 1 + 59 2 ) (A18) 

Here, /3i 2 = i^i 2 9 2 — [i 2 \9\. The corresponding expressions for r 2 i, z/ 22 , and v 2 \ can be inferred from Eqs. I|A12|) - (|A18|) 
by interchanging 1 *-> 2. For elastic collisions, the expressions (1A12|) — ( |A18|) reduce to those obtained for hard sphere 



The solution to Eq. (|A6|) is 

X a = (A- l ) aa ,Y (T ,. (A19) 

This relation provides an explicit expression for the coefficients <i", £i, and A^ in terms of the coefficients of restitution 
and the parameters of the mixture. Their explicit forms are 

d'l = 1 {2 [2v l2 Y 2 - Y x (2v 22 - 3C*)] [vnvzi - vnv 22 + 2(m + v 22 )C - 4C* 2 ] 



-2 (JQ (Y 3 + Y 6 ) [2v 12 v 21 + 2v 2 22 - Q*{7v 22 - 6C*)] 

-21/12 (f^l (n + ^e) (2^11 + 2^ 22 -7C)i, (A20) 



9xi 



h = {-2Y 3 [2{v 12 v 21 - ^11^22)^22 + C,*{7vuv 22 - f>v\ 2 v 21 + 2v% 2 - 6z/n£* - 7v 22 C + 6C* 2 )] 
+2Y A v 12 [2v 12 v 21 - 2v lx v 22 + 2C*(^n + ^ 22 ) - C* 2 ] 

+2Y 5 (* [2v 12 v 21 + v 22 (2v 22 - 7C) + 6C* 2 ] - 2is 12 CY 6 [2(v n + v 22 ) - 7C*]} , (A21) 

Ai = ^ {-Y 3 C [2v l2V21 + v 22 {2v 22 - 7C) + 6C* 2 ] + vuCY* [2(m + ^22) - 7C*] 

-Vs [4^i 2 ^ 2 i(^ 22 - C*) + 2^ 2 2 (5C* ~ 2i/ n ) + 2^ u (7i/22C* - 6C* 2 ) + 5C* 2 (6C* - 7^ 22 )] 

+1/12*6 [4^21 + 2^i(5C* - 2^22) + C*(10^22 - 23C*)]} , (A22) 



where 



A = [4(i/ 12 z/ 2 i - i/ii^ 2 ) + 6C*(Mi + K22) - 9C* 2 ] [^12^21 - ^22 + 2C>ii + ^22) - 4C* 2 ] . (A23) 

The expressions for d 2 , £ 2 , and A2 can be obtained from Eqs. (|A20|) - (|A22ll by setting 1 <-> 2. From the above 
expressions one can easily get the transport coefficients D' , L and A from Eqs. Q22 [l -H24| ) . respectively They are 
functions of X\ but independent of temperature and pressure. 



The pressure tensor Pj~,i is given by 
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where r\ is the shear viscosity coefficient. Its expression can be written asi£ 

nT 

V = — {xiliVi + X2llm) , 



with 



272(2A 22 - C*) - 471 A 



12 



7i72 [C* 2 - 2C*(An + A 22 ) + 4(AnA 22 - A 12 A 21 )] ' 



(A24) 



(A25) 



(A26) 



1)2 



2 7l (2Aii-C*)-47 2 A 2 



7i72 [C* 2 - 2C*(A U + A 22 ) + 4(A n A 22 - A 12 A 21 )] ' 
The dimensionless quantities Ay are given b y 12 ' 13 

2 



An = 



—xi ( — 

8 



-1/2 



1-^(1 -an) 2 



+ yjw2i(l + a 12 )el /2 6 2 1/2 [60r*(/xi 2 2 - M2i^i)(^i + fe) 

+^2i^r 2 (^i + ^) 1/2 (3 - a 12 ) + s^ 1 ^ + 2 r 1/2 



-1/2 



(A27) 



(A28) 



A12 — 



8 „. M21/-, , „ s fl 3/2 fl -l/2 

15 /U.12 



60 2 - 2 (Mi2^ 2 - mhWi+h) 



-1/2 



+ 2^2i^ 2 (#i + ^ 2 ) 1/2 (3 - a 12 ) - 56^(01 + e 2 )-^ 



(A29) 



The corresponding expressions for A 22 and A 2 i can be inferred from Eqs. i|A26() and (|A27I) by interchanging 1 <-> 2. 

The program for calculating the cooling rates, the temperature ratio and the transport coefficients of the binary 
mixture can be obtained on request from the authors. 
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